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Abstract — The intermediate transport regime in nanoscale 
transistors between the fully ballistic case and the quasi- 
equilibrium case described by the drift-diffusion model is still an 
open modeling issue. Analytical approaches to the problem have 
been proposed, based on the introduction of a backscattering 
coefficient, or numerical approaches consisting in the Monte 
Carlo solution of the Boltzmann transport equation or in the 
introduction of dissipation in quantum transport descriptions. 
In this paper we propose a very simple analytical model to 
seamlessly cover the whole range of transport regimes in generic 
quasi-one dimensional field-effect transistors, and apply it to 
silicon nanowire transistors. The model is based on describing a 
generic transistor as a chain of ballistic nanowire transistors in 
series, or as the series of a ballistic transistor and a drift-diffusion 
transistor operating in the triode region. As an additional result, 
we find a relation between the mobility and the mean free path, 
that has deep consequences on the understanding of transport 
in nanoscale devices. 

Index Terms — nanowire transistors, quantum wires, ore- 
dimensional transistors, ballistic transport, compact model, drift- 
diffusion transport 



I. INTRODUCTION 

Multiple gate architectures such as gate-all-around (GAA) 
MOSFETs have lately attracted significant interest [1], [2], [3], 
and have emerged as promising options to keep short chan- 
nel effects under control, exhibiting quasi-ideal subthreshold 
swing with undoped channels. This has the very important 
consequence of alleviating intrinsic variability of transistor 
threshold voltage, which in planar MOSFETs is mainly due 
to channel doping. 

Nanowire FETs are a particular case of multiple gate FETs, 
in which quantum confinement occurs in the transverse cross 
section of only few nanometers. Nanowire FETs are basically 
quasi one-dimensional transistors, where transport occurs in a 
set of loosely coupled propagating modes. 

From the point of view of modeling, several papers have 
appeared in the literature addressing transport and quantum 
confinement in silicon nanowire transistors. In pioneering 
works [4], [5], [6], [7], the effects of quantum confinement 
on a silicon nanowire were discussed. Numerical detailed 
investigations of quantum confinement in silicon and silicon 
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germanium nanowires with the anisotropic effective mass 
approximation, and its effect in lifting the degeneracy of 
silicon conduction band minima were discussed in [8], [9]. The 
electrostatics of silicon nanowire devices with cylindrical sym- 
metry has been investigated through a perturbative approach to 
the Schrodinger equation [10] or with a self-consistent solution 
in Poisson-Schrodinger equation with cylindrical coordinates 
[11]. 

Analytical models of ballistic nanowire transistors have 
been proposed in [12], [13], [14] and a broad review can 
be found in [15]. In real nanowire devices currents are much 
lower than those predicted by ballistic models [2], which can 
only be used as an asymptotic performance limit. 

Non ballistic transport in quasi- ID channels is harder to 
model. As far as numerical studies are concerned, far-from- 
equilibrium transport in silicon nanowire transistors was in- 
vestigated in [16] within the non-equilibrium Green's func- 
tions formalism, for both ballistic and dissipative transport, 
using the Biittiker probes approach to model inelastic scat- 
tering. A subband-based drift-diffusion simulation, in which 
the 3D electrostatics is solved self-consistently with the 2D 
Schrodinger equation in each transverse cross section and a 
set of ID continuity equations based on the drift-diffusion 
description, has been proposed in [17]. 

As far as analytical models of dissipative transport in quasi- 
1D FETs are concerned, notable examples are Ref. [13], which 
proposed a semiclassical model with drift-diffusion transport 
and constant mobility inside a cylindrical MOSFET, and Ref. 
[18], in which a polynomial expansion of the Fermi integrals 
for the mobile charge is used. 

Specific scattering mechanisms such as phonon scattering 
have been numerically addressed within the non-equilibrium 
Green's functions approach by Jin et al. [19] and by M. Gilbert 
et al. [20], [21]. 

We believe it would be very interesting to have an analytical 
model capable to seamlessly cover the continuum of transport 
regimes between the limits of ballistic transport and drift- 
diffusion (i.e. quasi-equilibrium) transport. Such a model, 
theoretically derived from the formalism of Biittiker virtual 
probes [22], and consisting in either a chain of ballistic 
transistors or in the series of a fully ballistic and an ideal 
drift-diffusion transistor, has been proposed in [23], [24] for 
2D MOSFETs. 

On the basis of this work, in the present paper we present 
an analytical model capable of describing the complete range 
of transport regimes in quasi- ID FETs, from fully ballistic 
to long channel quasi-equilibrium drift-diffusion behavior. A 
preliminary attempt has been proposed in [25]. As we shall 
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Fig. 1, Schematic band diagram of a ID SNWT. V g is the gate Fermi 
potential, 0,„ is the gate workfunction, x is the channel electron affinity, C 
is the potential in the centroid layer. qe a are the eigenvalues of the vertical 
confinement, drawn out of scale. 

show, the model is sufficiently simple to be suitable for circuit- 
level simulations and provides a strong intuitive picture of 
the transition from ballistic to drift-diffusion transport, which 
is missing in other descriptions of partially ballistic transport 
such as those relying on the introduction of a backscattering 
coefficient [26], [27]. 

The paper is organized as follows: in Section II we set up 
a model for ballistic transport in a nanowire transistor, that, 
in Section EH, we apply to the case of a chain of ballistic 
transistors. With a linearization procedure we show that a 
sufficiently long chain of ballistic transistors, can be regarded 
as a drift-diffusion channel. But the same approach fails for 
short ballistic chains, in which transport has an intermediate 
nature between ballistic and drift-diffusion. This difficulty is 
tackled in Section IV, where we present a compact model 
for intermediate transport that treats the ballistic chain as 
a series of one drift-diffusion section, for the first N - I 
transistors, and the remaining one ballistic channel, in which 
the non-equilibrium character of the intermediate transport 
manifests itself. In Section V we introduce the cylindrical and 
rectangular confinements for silicon nanowires considered in 
this paper and in Section VI we compare the results of our 
drift-diffusion and intermediate transport compact models with 
the numerical solution of the transport through ballistic chains 
of different length. In the section we also give an estimation 
of the current ballisticity ratio as a function of the transistor 
chain length. 

H. BALLISTIC TRANSPORT 

In the following discussion, we describe our approach in 
the general situation of a n-FET with a quasi- ID channel, 
with the effective mass approximation. Indeed the subband 
energies are determined by the transverse confinement, and 
to explicitly account for the capacitive coupling between gate 
and channel, the contact geometry has to be taken in account. 
The bottom of the ID conduction subbands are formally 
defined as the sum qe a - q<p c , of the eigenstates of the vertical 
confinement with respect to the conduction band edge in the 
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Fig. 2. Schematic capacitance diagram for a wire with square and circular 
cross section, representing the series of the oxide capacitance C ox and the 
silicon capacitance Cj. In the square shape the interface between the silicon 
and the insulator is considered approximately isopotential and the mobile 
charge layer (the centroid zi) is approximated with a square contour. The 
table summarizes the geometrical parameters used in the simulation of the 
rectangular wire (R) and cylindrical one (C). 

centroid layer (qs a ) and of the electrostatic potential energy 
(-q4> c ) in the centroid layer, that is where one can think all 
charge localized, following the approach of Refs. [28], [23]. 
For simplicity, a denotes the set of the quantum numbers 
specifying the confinement. The dimensionality also modifies 
the Fermi-Dirac integrals F v -i/2 and F v entering the ballistic 
equations for the mobile charge and the current in the channel, 
respectively. In particular for a ID conductor, in effective mass 
approximation, v = 0, whereas for a 2D MOSFET v = 1/2 
[24]. A definition for the Fermi integrals, with v > -1, is 

FM = r , 1 7T f * dx (1) 
r(v + 1) J e x n + 1 

o 

with T being the Gamma function, acting as a normalizing 
factor for the Fermi integrals. For v < — 1 we can rely to their 
property (d / 'drf)F \,(rj) = F v -\(rj) for their definition [29]. 

We start from a generic multi-subband degenerate version 
of the ballistic model in [24]. The vertical electrostatic model 
we propose, is similar to that in [12], [13], it is also somewhat 
less sophisticated, because we will suppose that screening can 
be considered constant. This is done in view of obtaining an 
analytical model of intermediate transport. Anyway we note 
that this assumption is sound enough for small cross-section 
channels and low electron densities [30]. Indeed, in the inverse 
layer centroid approach, we consider the charge accumulated 
in the centroid layer and its geometrical screening is included 
in the effective gate capacitance as a series contribution Cd, 
therefore the effective gate oxide capacitance for unit length 
is given by: 



as reported in Figure [2] where the expression C ox depends on 
the geometry. 

If we suppose an undoped channel, consistently with FigQ] 
the linear density of mobile charge on the peak of the potential 
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barrier in the channel is given by: 

Q m = -c s [v g - {fa -x)lq- <fic] , (3) 

where fa is the electrostatic potential in the centroid layer, 
(0m - X) /<? is me f at band potential, given by the difference 
between the gate workfunction fa and the channel electron 
affinity x- 

In the case of ballistic transport, there is no local equi- 
librium so that no quasi-Fermi level can be locally defined, 
because two different carrier populations exist, originating 
from source and drain, that can be considered at equilibrium 
with the injecting electrodes, as discussed in [31]. These two 
populations are separated by the peak of the barrier in the 
channel, that controls transport. Therefore, only three points 
are important: source, drain and the peak of the electrostatic 
potential. In ballistic transport, carriers move without inelastic 
scattering along the channel and therefore at the subband peak 
the carriers that propagate toward the drain ("forward states") 
only come from the source, whereas the carriers propagating 
toward the source ("reverse" states) come from the drain. As 
a consequence, we have the superposition of two hemi-Fermi- 
Dirac distributions. Following the considerations in [12], we 
can write for the ballistic mobile charge linear density on the 
peak: 

(4) 



(a) 



Qm = -qY^NalF-xninl) + ^-1/2(7$)], 



where: 



k B Tm a 1 



is one half of the effective density of states of the a-th 
subbands multiplied for its degeneration index g a . F^x/iiv) 
is the Fermi integral of order -1/2 and 



m = 
n a d = 



- V s - E a )/<f>t, 
~ Vd ~ S a )/fa. 



(5) 



The Fermi potential at the source (drain) is V s (j) and <p, = 
K/,T lq is the thermal potential. Equations 0) and © have to 
be solved simultaneously to obtain fa and Q m . 

From the Landauer formula, we can write the current as 
[12]: 

Id, = qY.G a [F (Tti)-Fo(T£!\, (6) 



where G a = gJ^£-T{\) is the effective injection rate of a ID 
channel multiplied by the degeneration g a of the a-th subband. 
Note that the current is dependent on the channel potential 
through r] s and rjd- 

III. FROM BALLISTIC TO DRIFT-DIFFUSION 
TRANSPORT 

We follow the approach developed in [23], [24] for a 2D 
MOSFET for the non-degenerate and degenerate cases. Here 
we analyze the case of a silicon nanowire transistor (SNWT) 
where the different dimensionality leads to different Fermi 
integrals entering the current and the charge expressions, and 
to different electrostatics. Moreover we considered a multi- 
subband degenerate model while in [24] only a single subband 
was presented. We obtain also a correcting factor for the 
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Fig. 3. (a) Circuit model of a generic SNWT, subject to inelastic scattering, 
in terms of a convenient chain of ballistic (B) SNWTs. (b) Approximate 
aggregation of the first N - 1 ballistic transistors in an equivalent Drift- 
Diffusion (DD) one. The macromodel DD+B comes out to be a suitable model 
for a device in intermediate transport regime. 



degenerate case, correcting the result of the linearization 
process whenever the low field approximation is not in fully 
satisfied. 

We recall that, within the Biittiker probes approach, inelastic 
scattering is thought as localized in special points, spaced by 
a length defined as "mean-free path" A. The virtual probes 
act as localized reservoirs along the channel, in which carriers 
are thermalized in equilibrium with their quasi-Fermi potential 
Vjfc, while transport from one virtual probe to the next is 
considered purely ballistic. We have a drift-diffusion transistor 
when the channel length is much longer than the free mean 
path, that from our point of view it is equivalent to have a 
long enough chain of ballistic transistors, as rigorously shown 
in [23]. On the contrary, when the number of internal contacts 
is small, transport is far-from-equilibrium, and is fully ballistic 
in the limit N = 1. We remark that within the Biittiker probe 
approach, transport of hot electrons is accounted only inside 
each ballistic channel, whereas a full thermalization occurs 
in correspondence of each probe, where electron density is 
described by a single quasi-Fermi level. It would also be 
interesting, but out of the scope of the present paper, to 
couple the transport equation with a heat diffusion equation, 
accounting for the energy losses in the Biittiker probes, leading 
to a non uniform temperature distribution in the device. 

We define Vj as the quasi-Fermi potential of the k-th virtual 
probe, and suppose that the k-th contact is placed at = kA 
with k = l,...,N, where the boundaries are fixed as Vo = 
V s = V and Vn = Vd — Vds- That is equivalent to place N 
ballistic SNWTs of channel length A in series, as sketched in 
Figure [3] Since the current h in any k — l,..,,N FET must 
be equal to Ids, we have N equations determining the local 
Fermi levels: 



h = qYG a [F {Tf k _ l )-F {nl)} 



(7) 



We note that rf k = (fa — — s a )/(p t where fa is the 
electrostatic self consistent potential in the conduction band 
peak of channel k, between the source contact k — 1 and the 
drain contact k. Introducing the definition V* = (Vk-i + V)t)/2, 
i.e. the mean potential between the two contacts (k and k— 1) 
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of channel k, and making explicit the Fermi integrals, we can 
rearrange (Q as follows 



CO 



-dx (8) 



cosh(^)] + [cosh(^)-l] 

where fjfc = k - V* - £ ff )M- 

At this point, in order to build an analytical model, we 
consider a large number of contacts and, having in mind 
that the current is constant along the channel, we extend 
Vk to a continuous quasi-Fermi potential V(x) satisfying the 
conditions: 



V 



^ X k +X k -l ^ _ 



dV ( x k + x k -i ^ _ V k - V k - 



(9) 

, ■ (10) 

ax \ 2 ) A 

Under the particular hypothesis that every ballistic SNWT 
works in the linear region i.e. that: 



dV(x) 

A—^- = V k - Vjm « 2<p„ 
dx 



(11) 



and expanding the terms sinh and cosh to the first order, we 
can put ([8} in the local form 



qA dV(x) v 

T a 

where the quantities pertain to the point x, as 

r, a (x) = [cf>(x)-V(x)-s a ycf >t . 



(12) 



(13) 



Evidently, if a voltage V ds is applied to the chain, ( fTTT i is 
satisfied if Vd s « 2<p,N. Moreover, within the same approxi- 
mations, the vertical electrostatics becomes: 

fi(*) = J] Q" {x) = 2q^N a F^ /2 [^(x)]. (14) 

a a 

An important aspect of ( flZb and (fl4l is that the current / a 
in subband a can be written in terms of the mobile charge 
density Q a (x) and of a mobility yf de (x) as: 

dV(jc) 



I"(x) = fi a deg (x)Q a (x)- 



dx 



(15) 



where the conduction is affected by the ID electron gas 
degeneracy through the mobility fideg(x). The mobility is given 
by 

„ v a A F-drfjx)] 

VdJ x ) = Tr~n r <u w ^ 

aeg 2<p t F-i/ 2 [lf(x)] 

where v a A/2(f>,, to which the expression ([Tol l is reduced in the 
non-degenerate limit, represents the low-field mobility for a 
ID gas of incoming electrons described by an hemi Maxwell- 
Boltzmann statistics [24] occupying the a-th subband, whose 
mean electron velocity is v a = ■jjjj^t characterized by a 
ballistic motion for paths of length / < A and a sudden 
and complete scattering at / = A. Considering the mean free 
path A as a constant, ( TTo*l > describes the degradation of carrier 
mobility due to degenerate conditions of Fermi-Dirac statistics. 
An analogous expression was recognized in a Monte Carlo 
simulation [32] for the case of strained silicon FETs. While, 



in [16] and [27], a similar, but not identical relation between 
the mean free path and the effective mobility has been found. 

The current is expressed in a local form in (fT2l and we can 
eliminate the gradient of the local quasi-Fermi level integrating 
along the channel and exploiting current continuity, leading to 



dV(x) 

(x)] —Ll dx . ( 17) 
dx 



In order to obtain a more compact form of ( fT71 ), we can change 
the integral variable as follows 

L V d rid 

r dv(x) r r dv 

qF. 1 [r,(x)]—^dx= F. 1 [v(V)]dV= F.^—dr], 
J dx J J dr/ 

V s 17, 

(18) 

where the term dV/drj is obtained by differentiating ( fT4b in dV 
and using the fact that drj/dV = (d(p/dV- l)/<p r . The current 
can be obtained with a numerical integration of 



ja 



qAG. 



,L J 



F-dV a (V)]dV, 



(19) 



where we note that r\ a not only explicitly depends on V, but 
also implicitly through <p c , as shown in ( fT3l ). Such dependence 
has to be taken into account in the self-consistent solution of 
the vertical electrostatics. 

Finally, we obtain the compact expression (which we will 
refer as the DD model) for the source-drain current for each 
subband a 



ZpPfiF-i.-wW, rf s ) - r-u- 3/2 (rf d , rf d )]). 



(20) 



Where for simplicity we defined 



r-i,-3/2(rf,Tf) = j ' F- X (x)F- V1 \x+ 8 -^^\dx, (21) 



with p a = ^j. It is worth noting that, as observed also in 
[23], [24], in the non-degenerate limit the first term in (120b 
reduces to the diffusion term of the EKV model, while in 
the degenerate limit, it corresponds to the current of a single 
ballistic channel divided by N. The second term of (f20b is 
instead associated with the drift current. 

We note that actually, in the low field approximation, for 
the case of a ID channel, the integral (O can be analytically 
solved as discussed in Appendix U However, the use of 
the analytical expression leads to the more complex and 
numerically expensive expression for the drift-diffusion current 
(|29] i, while giving only a slight correction of f!9[ . The model 
that employs the analytical solution of (|8), with the use of 
( |2"9i >, will be referred as DD* model. We conclude stressing 
the fact that DD and DD* models, considered alone, are not 
appropriate to describe transport whenever the condition (ITTb 
is not satisfied: for example in very short channels, where 
intermediate transport is expected. 




Fig. 4. Discrete quasi-Fermi potential for a chain of N = 5 SNWTs at fixed 
V,ij = 0.5 V. Fermi level is defined only on the virtual probes placed at points 
x = kX with k = 1 ... N. It is evident that the first N - 1 transistors work 
approximately in linear regime, while increasing V g a non-linear behavior is 
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Fig. 5. Comparison between NANOTCAD 2D simulation and the compact 
model calculation of the vertical electrostatic for the R and C SNWTs. The 
centroid potential obtained with NANOTCAD 2D (solid line) and the centroid 
potential obtained with the compact model (dotted) are reported. The compact 
model well accounts for the centroid potential, point in which the charge of 
the channel can thought to be concentrated. 



IV. COMPACT MODEL FOR INTERMEDIATE 
TRANSPORT 

Now we are interested in the development of a model that 
will be effective in the whole range of transport regimes, as 
proposed in [23] for the 2D MOSFET case. It is evident that 
in the general case of intermediate transport, the simplifying 
hypothesis (fTTT i. that enforces each SNWT of the ballistic 
chain to operate in the linear region, does not hold, and 
we can expect that some elementary channels can work in 
the saturation regime, or near it [23]. The behavior of a 
SNWT operating in such intermediate transport regime, can 
be obtained by solving the system for the complete ballistic 
chain 0, but it can represent a heavy computational burden, 
especially for a large number of internal nodes. In order to 




Fig. 6. Output characteristics for a chain of N ballistic cylindrical C SNWTs 
(denoted by NB), for N = 10, 20, 30, compared withe results of the DD 
model with A = L/N. The gate potential value is V g = 0.8 V. As the number 
of ballistic transistors in series increases the drift-diffusion approximation 
become more and more able to capture the behavior of the device. 



build a simple model, that can be more easily handled, we note 
in Fig. |4] where we plotted the quasi Fermi potential on the 
virtual probes for a chain decomposition of a SNWT, that when 
the saturating behavior of the elementary ballistic transistor 
emerges, it is present mainly on the last ballistic transistor of 
the chain. This non-linear behavior is a general condition for 
transistors in intermediate transport regime, due to the fact that 
in its end the channel narrows down and therefore, to maintain 
constant the current flux along it, a major spacing between the 
last source and drain levels is required (or, in the continuous 
limit, a steep drop of the quasi Fermi level). This fact suggests 
that we can aggregate the first N - I ballistic transistors in 
an approximate equivalent drift-diffusion transistor with ratio 
L/A = N—\ working in low field conditions, as it is represented 
in Fig. |3] and similarly to what was proposed in [23] for 
MOSFETs. We can therefore see a SNWT in intermediate 
regime as the series of a drift-diffusion (DD) channel for the 
first N - 1 transistor, that we can solve using ( fl9] >, and a 
single ballistic transistor, governed by (O. We only need to 
solve the DD+B system, imposing constant current through 
the two sections in series and solving the Fermi potential on 
the internal node between the DD and the B channels. We 
point out that in the DD+B model the ratio L/A has no need 
to be integer, because we apply the continuous DD equation 
(Ell. 



V. APPLICATION TO GAA-SILICON NANOWIRE 

At first we will consider a silicon wire with square cross 
section of side f „. Concerning the effective gate capacitance, 
unfortunately no simple analytical closed form is available. 
In order to simplify the electrostatics, we suppose that the 
interface between the oxide and silicon is approximately 
isopotential and we consider the variational approximation 
reported in [33], where the capacitance per unit length of a 



rectangular coaxial line, associated with the oxide layer, is 

8e ox 



Cox — 



ln(l + ^) 



(22) 



and similarly the capacitance associated with the charge in the 
silicon body C4 is 

8ns ci 

C d = — (23) 



The use of Cd, with Zi adequately chosen, permits to treat the 
capacitance due to the charge distribution in the channel cross 
section, in series with the oxide capacitance C ox . The term 
Zi is a characteristic dimension of the closed line where we 
can effectively localize the whole mobile charge. Here it is 
used as a fitting parameter for simplicity, while it is actually 
dependent, due to volume inversion, on the charge density in 
the channel. We point out though that it is smoothly varying 
for small section nanowires and low electron density [30]. 

In case of rectangular quantum confinement, the eigenvalues 
of the Schrodinger equation can be considered for simplicity 
[12]: 



ml?. 



(24) 



where the mass tensor can be defined as: 
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mi and m t are the longitudinal and transverse components ol 
the effective mass tensor of the degenerate minima of the 
conduction bands in Si. We can write the effective mass, 
characterizing the motion in the unconfined direction (z), as 
m a = m\ for a (100) silicon wire, with v running on the 
different Si conduction band minima. 

In the case of cylindrical quantum confinement, for the gate 
capacitance, we have instead: 

2nS °* (25) 



and 



r - 

^ox 



c d 



h(l + £) 



(26) 



Considering cylindrical quantum confinement [12], we have 
that the subband separation from the bottom of the conduction 
band is described by the approximated and handy expression: 



2„2 



2^ 



lR 2 



"1 + N - X - 



(27) 



where n\ is the radial quantum number, «2 the azimuthal 
quantum number and v runs on the different silicon valleys. 

We applied our model to a cylindrical quantum wire and 
a rectangular SNWT, denominated C, R respectively. Their 
geometries are described by Fig. [2] The inversion centroid 
layer depth was fixed in comparison with the 2D Schrodinger- 
Poisson simulator NANOTCAD [34], as shown in Fig. |5] 
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Fig. 7. Output characteristics of a C SNWT modeled as a chain of N = 5 
elementary transistors: 5B corresponds to the the exact numerical evaluation. 
The DD+B compact model is obtained considering the series of a DD channel 
governed by 120) plus a ballistic one and the DD*+B is analogous but uses 
)29t for the DD section. The choices of the gate potential are V g = 0.2, 0.4, 




Fig. 8. Output characteristics of the R SNWT modeled as a chain of N = 
2 elementary transistors: exact numerical evaluation (2B), compact model 
DD+B with a drift-diffusion plus a ballistic channels and DD*+B analogous 
to the latter except for the use of i29i in the DD section. The choices of the 
gate potential are V g = 0.2, 0.4, 0.6, 0.8 V. 



We note that a careful choice of Zi permits to recover the 
inverse layer centroid potential in full agreement with the 
NANOTCAD simulator, correctly accounting, thus, for the 
screening due to the charge inside the channel as a function 
of the gate potential. 

VI. RESULTS 

For a short channel transistor with length of few mean free 
paths, transport is quasi ballistic, and we have seen that the 
DD model ( fl9] > fails to describe its behavior. On the other 
hand, it is well known that the transport regime of a transistor 
with channel length much longer than the free mean path is 
described by the drift-diffusion model. We want to check if 
our model is able to correctly reproduce such transition and to 
investigate the number of free mean paths after which transport 
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Fig. 9. Transcharacteristic curves for a chain of N ballistic C SNWTs, with 
N ranging from 1 to 10 and Vj s = 0.5 V. The exact numerical evaluation 
(AT3), and the results of the DD+B and DD*+B compact models are shown. 
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L/X = N 

Fig. 10. Ballisticity index of a 7VB chain as a function of TV for the R, C 
SNWTs, and a MOSFET (see text for details). The gate potential is V g = 0.8V 
and the drain-source potential is Vj, = 0.5V. 



can be definitely associated to the drift-diffusion regime. In 
Fig|6] we plot the output characteristics for a chain of N 
ballistic channels numerically calculated (denoted Mi) and 
with the Drift-Diffusion approximation (DD characteristics), 
for N = 10, 20, 30 and gate potential V g = 0.8 V. We note 
that for a SNWT of length smaller than 10 A, a DD description 
is not appropriate. With increasing N the difference between 
the NB and DD models is reduced, and for a channel of length 
> 20A the output characteristics calculated with the DD model 
fully reproduce the corresponding numerically evaluated Mis. 

We have calculated the output characteristics of SWNTs 
described in the latter section employing a direct numerical 
solution of the chain of N ballistic transistors (NB), and 
compared them with our models for intermediate transport 
DD+B and DD*+B. Figure [7j and [8] respectively show the 
output characteristics for the C SWNT with N — L/A = 5, and 
for the rectangular SWNT with N - 2. Similar considerations 
apply to the two figures. While the DD approximated equation, 
derived from the linearization of the NB chain, inadequately 



reproduces the saturation behavior of the NB characteristics, 
the DD+B model seems suitable to describe SNWTs in the 
intermediate transport regime. As shown in UJ and [8] the DD+B 
and DD*+B models are really able to capture the non-linear 
behavior of the ATS transistors, although a non negligible error 
remains in the saturation regime. This is due to the the weakly 
non-linear transport in the DD section, that has been neglected. 
We note that in general the DD+B* improves the agreement 
with the ballistic chain characteristics. 

After testing our model with rectangular and cylindrical 
nanowire geometries, changing both the oxide and silicon 
length t ox , t s i and with different values of N, we can conclude 
that the DD+B and DD*+B compact model quite well repro- 
duces the output characteristics of degenerate SWNTs for any 
N, with errors in the saturation zone of few percentage points. 

In Fig. [9] the transfer characteristics of a C SNWT, treated 
as a chain of N elementary ballistic channels, with Vd s = 0.5 
V, is presented. Both the DD+B and the DD*+B models well 
reproduce the behavior of the corresponding ballistic chain in 
all gate voltage regimes, for all values of N. We note that 
the DD*+B model is always more accurate, in particular the 
correction is more evident for transistor with few nodes, at 
large gate voltage. 

We investigated also the so called ballisticity index of a 
transistor [35], that is given by the current ratio 7/4, between 
the current of the transistor and that of a corresponding 
ballistic one. The results of its calculation for a NB chain, 
with N ranging from 1 to 20, are shown in FigJTOl where 
we considered the R, C SNWTs and also, for comparison an 
undoped Double Gate MOSFET with t si = 4 nm, t ox = 2 
nm. The ballisticity index is monotonous and slowly decaying 
with N, the behavior is similar for all the transistors consid- 
ered here. The curves can be easily fitted with the function 
1/ [1 + r(N - 1)] where N is the number of ballistic elements 
and r ~ 0.25. We note that initially the ballisticity index 
steeply decreases with N. For longer channels the current 
becomes slowly varying with N: sliding from N = 10 to 
N - 20 the current only decreases from the 30% to the 20% 
of the ideal ballistic one. 

Having in mind a compact model, the calculation of equa- 
tions d2"0b or (fT9b for the DD+B model are still computation- 
ally expensive. Therefore we also tested the approximation 
of the integral in the drift-diffusion section (DD) with its 
symmetrical linearization [36], as discussed in Appendix ITT1 

VII. CONCLUSION 

We have presented a physics-based analytical model able 
to describe quasi one-dimensional field-effect transistors in 
the complete range of transport regimes extending from the 
fully ballistic case captured by the Natori model to the quasi- 
equilibrium case captured by the drift-diffusion description. 
Our proposed model sees a generic transistor as a long 
enough chain of elementary ballistic transistors in series with 
a common gate. Based on the Biittiker probes description of 
inelastic scattering, we have rigorously proved that the model 
reduces to the limit cases. In addition, as the most important 
result in this paper, we have shown that an equally adequate 



model, much simpler from the computational point of view, 
and more physically intuitive, is represented by the series oi 
an appropriate drift-diffusion one-dimensional transistor and 
a ballistic one-dimensional transistor, consistently with the 
results in [23], [24], that apply to 2DEG FETs. We have 
focused in this paper on silicon nanowire transistors, but om 
model is applicable without significant variations to any type 
of quasi-one dimensional FET, such as those based on carbon 
nanotubes, graphene nanoribbons, or other channel materials. 

Finally, we have shown that an interesting consequence oi 
our model is that, if a uniformly spaced chain is assumed, 
the Fermi-Dirac statistics degrades the low-field mobility, 
consistently with the observations in [32]. 
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Appendix I 

ANALYTICAL SOLUTION OF THE DD INTEGRAL 
We note that in the integral 



CO 

/: 
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in © has an analytical solution given by 
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where a = [cosh ^'J^' ) - lj /2. We replace the Fermi level 
difference between neighbor probes, by its mean value on the 
linearized chain 

V k -Vk-i * r = AV (DZ5) /(2<M0, 

where AV (DD) is the total potential drop in the DD section and 
N the number of elementary channels in it. 

In the low-field approximation we also replaced the 
sinh( ^20*"' ) w i m i ts arguments: we try to amend this by 
including a correction factor obtained by the ratio of the not- 
approximate term over approximate one 

'Vk-Vk-x\ ,Vk-Vk-i 



sinhf — „ , " ' I/- 



2<P, 



2<f>, 



(28) 



In the end, we reach a more accurate version of ( [T9| > for the 
DD section, given by the following expression 



r = 



qAG a sinh(y) 



tL yyja(a + 1) 
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+ 1) J 



I[V]dV 
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Appendix II 

SYMMETRICAL LINEARIZATION OF THE DD 
INTEGRAL IN SNWT 

The integral for the drift-diffusion current (T% is compu- 
tationally expensive for a compact model to be included in 
circuit simulators such as SPICE. Therefore we adopt a variant 



Fig. 11. Characteristic curves calculated for a R SNWT with N = 5. 
The results of the direct numerical calculation 5B, of the DD+B compact 
model and of the same compact model, employing symmetrical linearization 
DD/,„+B, discussed in the text, are shown. The curves are calculated for gate 
potential values of V g = 0.2, 0.4, 0.6, 0.8 V. 



of the symmetrical linearization [36], [37] in order to obtain 
an approximated result: 



J a 



<pc - V - s a dV X 

) dx 

d>, dx L 
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4>os 
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b c -V-e a ^dV X = 
cp, d(f> c L 

<f>, \d<Pc) m L 
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where we have defined the "quantum slope factor": 

1 



n q = 1 + 



4>c,m-V C j„-£„\ 



(30) 



that is a constant in the considered case. The linearization is 
done around: 

, <Pcs + <t>cd ,,,, 

<f>c,m = ^ ^ 3 ) 

where <f> cs and 4>cd can be obtained solving the vertical electro- 
statics dT4b . Moreover from vertical electrostatics we find V c m 
with an iterative process. We can observe in FigQl] that the 
symmetrical linearization of the DD integral well reproduce 
the not-approximated results. 
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